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Abstract: Fatigue describes the damage or failure of material under cyclic loading. Activation 
and deactivation operations of technical units are important examples in engineering where 
fatigue and especially low-cycle fatigue (LCF) play an essential role. A significant scatter in 
fatigue life for many materials results in the necessity of advanced probabilistic models for 
fatigue. Moreover, structural shape optimization is of increasing interest in engineering, where 
with respect to fatigue the cost functionals are motivated by their predictability for the integrity 
of the component after a certain number of load cycles. But mathematical properties such as 
the existence of the shape derivatives are desirable, too. Deterministic design philosophies 
that derive a predicted component life from the average life of the most loaded point on 
the component plus a safety factor accounting for the scatter band do not have this favorable 
property, as taking maxima is not a differentiable operation. In this work we present a new local 
probabilistic model for LCF. This model constitutes a new link between reliability statistics, 
shape optimization and structural analysis which considers the perspective of fatigue but also 
fits into the mathematical setting of shape optimization. The cost functionals derived in this 
way are too singular to be lower semi-continuous. We therefore have to modify the existence 
proof of optimal shapes [20] for the case of sufhciently smooth shapes using elliptic regularity, 
uniform Schauder estimates and compactness of certain subsets in C^{^'^^^ via the Arzela- 
Ascoli theorem, where Q^^^ is some shape containing all admissible shapes. Here, we analyze a 
class of cost functionals which consists of volume and surface integrals whose integrands include 
derivatives of the displacement field up to second order. Moreover, we extend our existence 
results to high-cycle fatigue (HCF) and deterministic models of fatigue. 



1 Introduction 



A design being made from material will fail if the material degradation due to loading exceeds 
certain limits. Reliability, i.e. the absence of failure, is thus the ultimate goal of structural 
design. Whether a design will operate safely under certain load conditions depends on the 
failure mechanisms which are very diverse for different material classes. Degradation can 
occur as a function of operating time, which is e.g. the case for creep damage. Or it can occur 
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1 INTRODUCTION 



when small plastic deformations under cyclic loading pile up and result in a crack. This is 
called fatigue which can be differentiated into high-cycle fatigue (HCF) and low-cycle fatigue 
(LCF). In this work we focus on LCF in conjunction with polycrystalline metals and consider 
the number of load cycles until crack initiation. 

The LCF failure mechanism - confer [B], [30], |15j and [HU] - with respect to polycrystalline 
metal leads to the stochastic nature of LCF which prevails on the macroscopic scale, where a 
statistical scatter of a factor 10 between the highest and lowest load cycles to crack initiation is 
a common phenomenon, even under lab conditions. This makes it obvious that the reliability 
of a design has to be understood quantitatively as the probability of the absence of failure 
after a given number of loads. Thus reliability statistics, see e.g. [13], should play a vital 
part in LCF design. Furthermore, LCF cracks initiate at the surface and they are small in 
the initiation phase. Therefore, they will only influence stress fields on micro- and mesoscales 
and so it is reasonable to assume that crack formation in one region of the component surface 
is not influenced by the crack forming process on another part. This forces us to see crack 
formation as a problem of spatial statistics [M]. We thus infer that hazard rates for crack 
initiation have to be integrals over some local function depending on the local stress or strain 
fields. Therefore, hazard rates for the component can be expressed by a surface integral over 
some crack formation intensity function depending on local fields. The latter can also be seen 
as the density for the intensity measure of a poisson point process (PPP) on [0, oo) x dQ, 
confer [25] and [5]. Here, the first component stands for 'time' measured in cycles. The very 
interesting article [15] also emphasizes the role of the PPP. 

In this work we model hazard rates and intensity densities of the PPP with a rather 
conservative approach. That is, we assume that the intensity measure is of scale-shape type, 
as it is the case for many distributions in reliability statistics. Furthermore, we assume that the 
scale variable Nf^^.^ is of the same functional form as the usual Coffin-Manson-Basquin equation, 
confer [6] and ([5]) below. Although we are not aware of any empirical study showing that 
LCF failure time distributions are Weibull, we determine the shape of the intensity measure 
according to this approach for mathematical simplicity. It then follows that the number N 
of cycles of first crack initiation will be a Weibull distribution as well. The Weibull scale 
parameter is given by the inverse of the L'"((90, IR)-norm of the inverse of variable iVjet- In 
contrast to [15] this purely phenomenological approach avoids detailed modeling at the meso 
scale which facilitates calibration with experiments. Both approaches have the use of Weibull 
distributions in common. From a materials engineering point of view our new model has the 
significant advantage - compared to standard methods in fatigue - of bypassing the standard 
specimen approach and of considering size effects. 

As LCF is driven by cyclic loads, the reaction of the component to these loads needs to be 
taken into account as well. Usually, this is done via a continuum mechanics approach and here 
we do follow. Furthermore, we restrict to linear elasticity as the fundamental equation that 
determines the loads. This seems to be in contradiction with the need to model also plastic 
deformation to capture the LCF damage mechanism. But there are time honored procedures 
that bypass this problem, namely elastic-plastic strain conversion of Neuber and Glinka, see 
[25] , Although those seems not quite up to date, they are still often used in engineering and 
they have considerable mathematical advantages when it comes to shape optimization - and 
this will be the next topic. 

So far we have considered the component shape Q as fixed. It is the task of design teams, 
to choose the shape such that reliability, performance and cost are in an optimal equilibrium. 
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Here, we consider the disciplinary task of structural design alone and hence focus on reliability. 
Narrowing in even further, reliability is understood as the probability of survival (no crack 
initiation) with respect to LCF at a given number of cycles. An optimal design with respect to 
LCF will then be a design, that maximizes, under given load conditions, the survival function 
SNin) = l — FN{n) with Fj^ the cumulative density function of the random LCF crack initiation 
time of the component. 

From a mathematical point of view, there are quite a few studies on the question of optimal 
design, see e.g. [20] and references therein. The problem of optimal design can be seen 
as a PDF restricted optimization problem with control parameter 0. Compactness of the 
set of admissible shapes and lower semi-continuity of cost functionals J{Q,u{il.)) to be 
minimized are essential ingredients. Here, n(il) uniquely solves the linear elasticity equation 
on with boundary conditions on dO,. In this context, linear elasticity is treated within the 
weak formulation [10]. Consider a sequence given by ti„ = u{^ln) G [H^{^ln)]^ for n G N. 
Under a few assumptions on dQn it has been shown that u'^^ — t- strongly in 
for a subsequence Un^, whenever ri„ Q in a suitable sense j20j . Here, contains all 
admissible shapes and u^^^ denotes a suitable extension of functions to il^^^ which are defined 
on a corresponding admissible domain. 

Unfortunately, in the context of LCF a problem occurs that prohibits the application of 
these results: The LCF failure mechanism leads to cost functionals J{il., u) - the probability of 
failure at given cycles with displacement field u - that depend on Vu and are surface integrals. 
While the components of Vn are in L^(il) and of u \qq in L^(9f2) for u £ [H^{fl)]'^, this does 
not guarantee that Vu is a well defined function. A proof of lower semi-continuity for the 
probability of crack initiation as a functional on therefore seems to be problematic 

or even wrong. Therefore, we have to look for an alternative strategy to prove the existence 
of shapes with optimal reliability for LCF. 

We address a general class of cost functionals which considers LCF survival and failure 
probabilities as defined in Definition l4. II which are apparently [C^(0''^*)]^ continuous. Although 
this seems to be a very rough way of modeling again, we impose smoothness conditions on dfl 
and on the loads and boundary values in order to obtain strong solutions in [C^(f2)]^ for a 
domain in the case of a fixed Dirichlet and a variable von Neumann boundary. Note that we 
even consider third derivatives to be able to extend our models to stress gradients. The key 
here is uniform elliptic regularity and Schauder estimates [2]. This is however still not enough 
as we need to prove the existence of a limiting function u £ [C^(r2°^*)]^ for a subsequence 
of solutions u{Qn) of the linear elasticity problem on admissible shapes suitably 
continued to Q^^^, Qn ^ But using the Arzela-Ascoli theorem, revisiting the Schauder 

estimates and proving that they are uniform with respect to the set of admissible shapes leads 
to the existence of such a limiting function. 

The existence of shapes that are optimally reliable with respect to the probabilistic LCF 
damage mechanism can thus be proven mathematically in the given set up and HCF failure 
mechanism can be considered as well. But this is not special for the probabilistic approach. If 
we apply a deterministic LCF lifing model, we need to minimize peak stress under given loads 
in order to maximise the LCF life. Peak stress is C^(Q^^^) continuous as well and thus we 
obtained existence of optimal shapes for the deterministic approach to LCF as a by-product. 
The same can be achieved for HCF. It is special to the probabilistic approach to fatigue that 
it is described by a cost functional which is a local integral of the displacement field and its 
derivatives. And this is exactly what is needed in shape sensitivity analysis j35j which we will 
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2 LINEAR ISOTROPIC ELASTICITY AND PLASTICITY 



address in our work |19j with respect to the probabilistic model. Shape sensitivity analysis 
in conjunction with the adjoint method is a highly effective procedure to find - after suitable 
discretization - shapes of (local) optimal reliability via a shape gradient decent method. Confer 
|32j for an application in a different context. 

This paper unifies some aspects of materials engineering, reliability statistics, elliptic PDE 
and shape optimization. We therefore give some compilations of known results from the re- 
spective fields to keep the article self contained. In particular, Sections [2] and [3] repeat some 
basics from materials engineering. Section H] motivates our main definition - the cost functional 
for local probabilistic LCF. Here, HCF in a probabilistic sense and a deterministic model for 
fatigue are discussed as well. An abstract setting for existence analysis of shape optimization is 
given in Section [5] where we also address C'^-admissible domains. This abstract setting provides 
the structure of Section [6] where we derive existence results for shape optimization problems 
with respect to linear elasticity and to C'^-admissible domains. Here, we focus on a general 
class of cost functionals of integral form without convexity constraints. At the end of Section 
[6] we apply our existence results to the probabilistic and deterministic models of Section H] for 
LCF and HCF. 

2 Linear Isotropic Elasticity and Plasticity 

In order to describe the behavior of components from poly crystalline metal under external 
loading, we employ linear isotropic elasticity which is justified by the assumption of isotropic 
material behavior at scales significantly larger than the grain size and of sufficiently small 
deformations. We closely follow [12], [23] and [6], respectively. 

Let C be a domain which represents the component shape filled with a deformable 
medium such as polycrystalline metal which is initially at equilibrium. Moreover, let be the 
normal of the boundary d^l, let / : — t- be an external load and let ti : — t- be the 
displacement field. Finally, let = dQo U di^N be a partition where dilu is clamped and on 
ddj\f a normal load g : d^l]\f — )• R'^ is imposed. Then, according to [12] the mixed problem of 
linear isotropic elasticity is described by: 

V -a^u) + f = inn, 

a'^iu) = X(y ■ u)I + fi{Vu + Vu^) inn, 

n = on dn^), 
a^{u) ■ V = g on SQtv- 

Here, A and /i are the Lame coefficients. The linearized strain rate tensor e^iu) : n — )• M^^^ 
is defined as e^(n) = ^(Vn + Vu^), i.e. efj = \ + for i,j = 1,2,3. Approximate 
numerical solutions can be computed by a finite element approach, confer [23] and |12j . 

In the following we phenomenologically discuss time-independent plasticity as we cannot 
assume a completely linear elastic material behavior when talking about fatigue. Plastic defor- 
mations can allude to an imminent residual fracture of the material. Therefore, knowledge of 
the threshold between elastic and plastic deformations is very important. According to Section 

3 in [^ this threshold is described by so-called yield criteria. For multi-axile stresses scalar 
comparison stresses ay are introduced. The comparison stress ay depends on the stress tensor 
(Tjj and points to beginning yielding if ay reaches a critical value cxcrit- So yield criteria can 
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be described by equations of the form y{(Jij) = 0, where y{(Tij) = (Jvi^ij) — <7crit- For negative 
values of y the material has solely elastic behavior. 

Plastic deformations within grains of a metallic material are caused by dislocation motions 
where crystal faces are displaced in different directions. Experiments have shown that only 
shear stresses can plastically deform metallic material. Thus, the deviation of the state of 
stress from the hydrostatic state of stress a' = a — ^tr{a) decides if metal begins to yield. 

In the space of principle stresses the yield surface circumscribes a cylinder around the 
hydrostatic axis cri = a2 = cr^. In this work we use the von Mises yield criterion which is given 
by a cylinder around the hydrostatic axis with a cross section of radius R = ^/2kp. Here, kp 
is the critical value of the von Mises yield criterion 



^tr(cT'2) = [(ai - C72)2 + {a, - a^)' + (^^2 - ^3)2] = kp. (2) 

The left-hand side is proportional to the elastic strain energy of distortion. If the criterion 
is applied to uniaxial tensile tests the relationship kp = Rp/^/S is obtained, where Rp is the 
critical value of the only nonzero principle stress in a uniaxial tensile test. Defining the von 

Mises stress as = ^tr{a"^) the criterion can be rewritten in the form = Rp and used 
to predict yielding of metal under any loading condition from results of uniaxial tensile tests. 

Now, we present the Ramberg-Osgood equation which is used to locally derive strain levels 
from scalar comparison stresses which determine strain-controlled fatigue life, confer |31] . This 
equation describes stress-strain curves of metals near their yield points and is very accurate 
in the case of smooth elastic-plastic transitions. The latter can be observed for metals that 
harden with plastic deformations, for example. Let K denote the strain hardening coefficient 
and n the strain hardening exponent. Then, the Ramberg-Osgood equation is given by 



E \K 

with Young' modulus E = ^^^^^^^^ and the equation defines the comparison strain £y. We 
also write = R0((T^). 

At the end of this section we introduce the method of stress shakedown by Neuber which 
describes a conversion of elastic to elastic-plastic stress, confer [28] and [6]. If linear elasticity 
leads to stress values which are greater than the material yield strength such a stress shakedown 
can be performed to compute elastic-plastic stress values. The stress shakedown is based on 
an energy-conservation ansatz which results in a relationship between the elastic von Mises 
stres^ af, and the elastic-plastic von Mises stress : 

- o-^, e^, = — + a^, — . (4) 



E " " E ' " \K 

Here, the Ramberg-Osgood approach is used. Kt is the notch factor, which is set to one if 
is obtained from the solution of a linear elastic equation ([T|) where notches are part of the 
boundary definition. Given the elastic comparison stress cr^, we can thus calculate the elastic- 
plastic von Mises stress by solving ^ and thus we are able to obtain Sy from ([3]). We also 
write ay = SD~^(o"^). 



■^Confer [21] for details on Neuber shakedown in conjunction with equivalent stresses. As an alternative to 
Neuber's rule we could have also used Glinka's method, confer e.g. [26j . 



6 



3 FATIGUE 



£a(log) 




Figure 1: Triangle shaped load-time-curve. Figure 2: EN-diagram of a standardized speci- 
men. 



3 Fatigue 



In structural analysis, fatigue describes the damage or failure of material under cyclic loading, 
confer [6] and [30]. Compared to the static case material is damaged by much lower load 
amplitudes of cyclic loading. Moreover, in ductile materials under cyclic loading there are less 
significant plastic deformations before failure so that an imminent damage is more difficult 
to detect. In many engineering cases the functional description of the loading cycles is very 
complex. This leads to a preference on constrained analyzes of representative cases of cyclic 
loading. Figure [1] shows a triangle shaped uniaxial load-time-curve as an example, where 
fa = (o"max " crmin)/2 is the stress amplitude. In material science standardized specimens 
under such uniaxial load-time-curves are analyzed to predict the behavior of more complex 
engineering parts. 

For backgrounds on surface driven LCF failure mechanism with respect to polycrystalline 
metal we refer to [B], [3D], [TS] and [5B]. Furthermore, recall our introduction. We now discuss 
models to describe failure times in fatigue. In fatigue specimen testing the number of cycles 
until failure is determined and if the tests are strain controlled so-called E — N diagrams are 
created, see the test points in Figure [2j The functional relationship of the corresponding strain 
amplitude and number of cycles is called Wohler curve. For the purpose of application and 
analysis of fatigue one subdivides the range of cycles in low-cycle fatigue (LCF) and high-cycle 
fatigue (HCF) although there is no sharp threshold between them, confer [33] . Note that LCF 
loads are often strain controlled. However, HCF loads are mainly stress controlled so that 
corresponding S — N diagrams are analyzed. 

Figure [2] shows the relationship of strain amplitude and the life time Ni to crack initiation 
measured in cycles. In the LCF range the plastic part has the greater contribution to the 
strain amplitude whereas in the HCF range the elastic part e^^ dominates. It is assumed that the 
additive law Sa = + £^ holds. The elastic part can be described with the Basquin equation 

= -^{2Ni)^, where a'j is the fatigue strength coefficient, b is the fatigue strength exponent 

and E is Young's modulus. For the LCF range the Coffin-Manson equation = £'f{2NiY 
describes the denominating plastic part, where £[ and c are empirical constants called fatigue 
ductility coefficient and exponent, respectively. A discussion of the physical origin of this 



7 



equation can be found in [36]. One finally obtains the Coffin-Manson-Basquin (CMB) equation 

ea = et + e^' = ^m)' + e'f{2N,r. (5) 

This equation can be calibrated according to the test points by means of maximum likelihood 
methods, for example. For this and other methods confer [T7| and |33| . 

Having analyzed standardized specimens one assigns that Wohler curve to an engineering 
part under cyclic loading which corresponds to the loading and temperature conditions at the 
part's surface position of highest stress. Safety factors are additionally imposed to account for 
the stochastic nature of LCF, size effect^ and for uncertainties in the stress and temperature 
fields. In the following we refer to this method as the standard method in fatigue which is 
widely used in engineering, confer [30], [B] and [55] . 

Note that several extensions exist to ([5]). In particular, when approaching the HCF region 
at ~ 10'' - 10^ cycles, mean stress effects are of increasing importance. The modified Morrow 
equation is one approach that would fit into the approach taken here without changing much 
of the mathematics. We refer to [30] for further discussions on that topic. 

4 From Reliability Statistics to Probabilistic LCF 

Since fatigue tests show a great scatter in the life time of specimens several fatigue tests have 
to be conducted to obtain statistical statements for life times. An overview over fatigue tests 
and statistical methods can be found in [33j and [3U]. In the following we employ the concept 
of reliability to derive our probabilistic LCF model. From a materials engineering and mathe- 
matical point of view our model has significant advantages compared to the standard method 
in fatigue, recall Section [3l On the one hand it bypasses the standard specimen approach and 
considers size effects, i.e. results from arbitrary geometries under LCF failure mechanism can 
be used to calibrate our model and every position of the surface of an engineering part is taken 
into account by a surface integral which does not need information on Wohler curves of a 
specific specimen. Thereby inhomogeneous stress fields are considered. On the other hand our 
model fits into the setting of sensitivity analysis in shape optimization what we will address 
in [12]. 

We model failure-time processes on continuous scale, although time in our context is a 
number of load cycles and strictly speaking is an integer number. Let denote a continuous 
random variable on some probability space with values in which describes the time of 
crack initiation here identified with failure of a system or component. Let P be the underlying 
probability measure and recall the cumulative distribution function F]\f{n) = P{N < n), the 
density function /iv(n) = dFj\f{n)/dn, the hazard function h{n) = /7v(n)/(l — F/v(n)), the 
cumulative hazard function H{n) = h[s)ds and the survival function S'Ar(n) = 1 — -F/v(n) = 
exp(— if(n)), confer [T^ . 

Let an LCF failure mechanism on a boundary of a domain C M'^ be given whose failure- 
time is described by the random variable N. Moreover, let {^i}i=i,...,m be an arbitrary partition 
of do,, i.e. dQ = L)"^iAi with Ai pairwise disjoint. Now, we introduce and discuss the crucial 
assumption for our local and probabilistic model of LCF in the case of polycrystalline metal: 



*Note that different geometries of test specimens lead to different Wohler curves, confer |30) . 
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4 FROM RELIABILITY STATISTICS TO PROBABILISTIC LCF 



Assumption (H) 

The LCF failure mechanism on di} induces failure mechanisms on each Ai with crack initiation 
times Ni for i = 1, . . . ,m such that the random variables {A^j}i=i,...,m are independent. 

The time of initiation of the first LCF surface crack on dO, is always N = min(A'^i, . . . , Nm)- 
This is identified as the cycles to failure for the component. If hi denotes the hazard function 
of Ni for i = 1, . . . ,m and assumption (H) holds it follows that h = "^^i hi is the hazard 
function of A^. That is why we call the previous assumption the property of spatial additivity 
of hazard rates. For polycrystalline metal assumption (H) will hold to a good approximation 
if the surface zone that is affected from the crack initiation process of a single LCF crack is 
small with respect to the surface of the component. This surface zone corresponds to faces of 
a few grains. As long range order phenomena are unusual in polycrystalline metal, we pass to 
the following stronger assumption: 

Assumption (L) 

In any measurable surface region A C 90 the corresponding hazard rate Ha is a local functional 
of the elastic displacement field u in that particular region, 

/i^(n) = / p{n;Vu,V'^u)dA. (6) 
Ja 

Here, Vu is the Jacobian matrix and V^u the Hessian of u. Note that the loads that we 
consider in the given context (mainly elastic plastic stresses and strains) can all be expressed 
as functions of Vu. 

Let us motivate that assumption (L) is only slightly stronger than assumption (H). Consider 
a surface portion A C over which the variation of loads is negligible. Then we can further 
subdivide this region into m smaller pieces Aj of equal surface volume \Aj\ = \A\/m. As the 
loads are approximately constant over the Aj, we should have h^j ~ hAi, hence the hazard 
functions for crack initiation in the respective portion of the surface are essentially equal. Since 
for whatever subdivision hA = ^JLi^^i assumption (H), we get that hA ~ rn x hAi = 

jjj-^hAi- For m — 7- oo we thus see that the limit p = limm-i-oo exists and one obtains 

h K, A - p. Here, p stands for a hazard density function. 

For mathematical simplicity, we restrict ourselves to most simplistic form of assumption 
(L), namely the dependence on elastic strains (or equivalently stresses), only. Modeling damage 
times in the spirit of elastic support factors depending on V^u via x* = l^iii might be an 
interesting option for the future. 

In the case of inhomogeneous strain fields assumption (L) implies the previous ansatz to the 
form h{n) = Jq^ p(n; e^) dA for some hazard density function p and where is the linearized 
strain rate tensor. Finally, we obtain the cumulative distribution function which yields the 
probability of failure in di} until cycle n: 



Fjv(n) = 1 - exp(-if(n)) = 1 - exp ^- ^ h{t) dt^ 

= 1 - exp ( - /" [ p{t; e"") dAdt] . 
\ Jo Jan J 



(7) 



The ansatz ([7]) can also be derived by the Poisson point process on [0, oo] x with p{n; e^{x)) 
as the intensity measure. The advantage of this point of view is that also the probability of 
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the occurrence of a given number of cracks initiated in ^4 C d^l within n load cycles can be 
calculated via the Poisson statistics as 

P(number of crack initiations on A = q) = e ^— for z = I I p{t;e'^) dAdt. (8) 

Q- Jo Ja 

But this approach will break down, if cracks have grown sufficiently large in order to mutually 
influence their local stress fields. For an introduction into point processes confer [3] and |25| . 

Not having specified the hazard density function p we now establish a link to deterministic 
LCF analysis which leads to an appropriate choice for p. Recall the CMB equation ([5]) which 
determines the functional dependency of strain and cycle of failure for a standard specimen of 
a fatigue test, confer Section [3] and [30]. The relation ([5]) refers to a specimen with a specific 
geometry in a homogeneous strain field. In order to convey our approach so far to an arbitrary 
geometry in inhomogeneous strain fields one has to interpret the CMB ansatz in a new way, 
but such that the results for standard specimens under homogeneous conditions are included. 

In order to be concrete, we will make the assumption that life cycles to crack initiation N 
are Weibull distributed: 

/ \ m—l 

p{n;x) = p{n;£%x)) = ^. .. I ,^ . I . (9) 



Here, A'dct(£'^) is the scale parameter and is computed as follows: Solving the elasticity problem 
([1]) results in the linearized strain rate tensor and the linear elastic von Mises stress field 
0"^. Applying the method of stress shakedown by Neuber to yields the elastic-plastic von 
Mises stress amplitude a a corresponding to uniaxial fatigue test. Using the Ramberg-Osgood 
equation yields the corresponding strain field Ea which is the left-hand side of the CMB equation 
([5]). The latter equation finally leads to N^{e^). The shape parameter m G (0, oo) controls 
the scatter with small m > 1 correspondqj to large scatter and the limit m — )• oo is the 
deterministic limit. The Weibull model is not necessarily true for real data. However, the 
Weibull hazard rate can be easily replaced by any other differentiable hazard rate with scale 
parameter A'det; without changing much of the content of this paper. Combining ([7]) and ^ 
yields 



F„(„) = 1 - exp I - £- ) ' aAi. I = 1 - exp 



n 



A^det 



(10) 



Definition 4.1 (Local and Probabilistic Model for LCF) 

Let C be a domain whose boundary dO. is subject to surface driven LCF failure mecha- 
nism. Moreover, let the scale field N^ct{x) = N(ict{.£a{x))i x E 90, be the solution of the CMB 
equation ([2]) 



Saix) = -^{2Ndct{x))'' + e'f{2Ndct{x)r, (11) 

where £a{x) is calculated from e'^{x) visH linear isotropic elasticity, from the von Mises stress 
ay{x), from aa{x) = SD~^{ai,{x)/2) according to the method of stress shakedown by Neuber 
and from the Ramberg-Osgood equation with £a{x) = RO(fTa(x)). Then, the local and proba- 
bilistic model for LCF is defined by the cumulative distribution function (jlOp for n G M.^ and 
some m > 1, which yields the probability for LCF crack initiation in the interval [0,n]. 



''O < m < 1 is not realistic for fatigue. 
^Confer Sections [2] and O 
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5 SHAPE OPTIMIZATION AND -ADMISSIBLE DOMAINS 



Note that CMB parameters in this approach are not the same as obtained from fitting 
standard specimen data. In the local, probabilistic approach to LCF the CMB parameters need 
to be fixed by the usual maximum likelihood methods of reliability statistics, confer [T7] and 
[13j . Volume driven HCF failure mechanism can be considered by replacing the surface integral 
in (jlOp with a volume integral whose integrand only differs by different material parameters. 
For a discussion of HCF failure mechanism confer [B] and [3D]. In case of mixed LCF and HCF 
failure mechanism a volume integral of similar form can added to (jlOp . 

In contrast to the probabilistic model the standard method in fatigue yields the following 
deterministic LCF model: 

Definition 4.2 (Deterministic Model for LCF) 

Given the setting of the previous Definition 14. II the deterministic model for LCF is determined 
by inf{A^det(x)|x G On}. 

In the next sections we will prove the existence of optimal shapes for Definitions 14.11 and 14.21 

5 Shape Optimization and C'^- Admissible Domains 

In the following we introduce into an abstract setting of shape optimization and discuss so- 
called C'^-admissible domains. This leads to a theoretical frame for an existence proof of 
optimal designs with respect to a general class of cost-functionals which include functionals 
of LCF and HCF failure mechanism in our probabilistic and deterministic sense. First, we 
address basic notations and closely follow Section 2.4 in [20]. For further introductions into 
shape optimization confer [35], [8] and [32], for example. 

Notation (Family of Admissible Domains, State Space) 
Let O denote a family of admissible domains and let V{0,) for every 17 € O denote a state 
space of real functions defined in $7. 

Notation (Convergence of Sets and of Functions with Variable Domains) 

Let {Qn)neN be a sequence in O and let Q, £ O. Then — > as n — t- oo denotes the 
convergence of {^ln)neN against Cl. If (yn)nGN is a sequence of functions with y.„ G V{Qn) for 
every n S N and if y G V{Cl) then y„ ^ y as n — )• oo denotes the convergence of (yn)nGN 
against y. Moreover, it is assumed that any subsequence of a convergent sequence converges 
against the limit of the original one. 

In every Q £ O one solves a state problem which can be a PDE or a variational inequality, 
for example. Assuming that every state problem has a unique solution and associating with 
any S O the corresponding unique solution u{Q) G V{Q) one obtains the map u : i — > 
u{n) e V{n). Let O be a subfamily of O, then G = {{n,u{n)) | G O} is called the graph of 
the mapping {u{-)) restricted to O. 

Definition 5.1 (Cost Functional, Optimal Shape Design Problem) 

A cost functional J on O is given by a map J : (f7,y) i — > J{Q,y) € M, where 17 G O and 
y E V{Cl). Let O be a subfamily of O and for every 17 E O let u{fl) be the unique solution of 
a state problem given in 17. An optimal shape design problem can then be defined by 



r Find 17* E O such that 

\ J(17*,u(17*)) < J(17,n(17)) V17 EC 



(12) 
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Now, we present a statement regarding the existence of optimal shapes. Note that this 
theorem will structure the following section where we prove of our main existence results. 

Theorem 5.2 (Existence of An Optimum in Shape Design Problems) 

Let O he a family of admissible domains and O a subfamily. Moreover, let J be a cost functional 
on O and assume that every Q G O has a state problem with state space V{p?) where each such 
state problem has a unique solution u{Q) G V{^). Finally, conjecture 

• Compactness of Q = | G O}: 

Every sequence {yt,n,u{yt,n))n^n C G has a subsequence (fi^^, u(r2nj.))„j,gN which satisfies 

u(J7„^) ^ u{Vt), A; — 7- 00 

for some (r2,n(il)) G Q. 

• Lower semi- continuity of J: 

Let {^ln)neN with r2„ G O, n G N, and (?/n)n6N with yn G V{Qn), n £ N, be sequences 



and let and y be some elements in O and in V{Q), respectively. Then 

=^ \mimin^ooJ{^n,yn)>J{^,y)- 



Q.n ^5 n — )■ 00, 
yn-^y, n 00 



Under these assumptions the optimal shape design problem il2\) possesses at least one solution. 

According to the previous theorem and abstract setting we have to define the family of 
admissible domains O. In this work we consider so-called C*^-admissible domains which have 
smooth boundaries. For the sake of simplicity we only optimize a part of the boundary of 
these shapes. This method is described in Section 2.8 of [35] and in [20], too. Importantly, 
C*^-admissible domains satisfy compactness properties as required in Theorem 15.21 according 
to the Arzela-Ascoli theorem. 

At first, admissible domains are defined which are determined by uniformly bounded func- 
tions. Later on, these functions are assumed to be sufficiently smooth to obtain C'^-admissible 
domains. On these domains we will later impose boundary value problems of linear elasticity 
which are characterized by disjoint Dirichlet and Neumann boundaries. In the following we 
employ so-called uniform cone properties which are defined in Appendix A of |20j . 

Definition 5.3 (Basic Design, Design Variables, Admissible Domains) 

Let C be a simply connected and bounded domain which is called basic design under the 
following assumptions: 

• Cl has a uniform cone property and a C'^-boundarjj^ for some A; G N, A; > 1. 

• There is some amin G M so that the cross section = {{xi,X2) G | {xi,X2, amin) G ^} 
is a nonempty domain in M^. 

^This leads to a local description of dCl by a finite number of hemisphere transformations of class C*^, confer 
Definitions I A. II and [X2I and Remark IX3l 
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5 SHAPE OPTIMIZATION AND C^^ -ADMISSIBLE DOMAINS 



3^3 




Figure 3: Admissible shape. 

• There are some z £ Cl and r > such that B{z, r) C 0, and Z3 + r < ttmin- 
For Omax > ckmin and positive constants Li, L2, L3 the elements of 



ad 



amiri < a < amaxin Om, a 



2di "|(9!^2 



ttmiri) / a(x)(ix = Li, ||a||(;;fc < ^2 



< Lsllx - y\\2 Vx,y e Osdj 

are called design variables. Let a G U"''^ define the set 

r2(a) = {x £ Cl\x3 < amin} \ B{z, r) U {x G | (xi, X2) G ^2^, amin < 2^3 < a(xi, X2)}, 

see FigureO Varying only functions a G IJ^ the domains O = {^{ot) \ a G IJ^^ with Lipschitz 
continuous boundaries are called admissible domains. Finally, choose a superset il'^^* such that 
extension Lemma I A. 41 can be applied. 

Lemma 5.4 IJ"-'^ is compact in [C^{Q2d)^ II ■ Wc'') ■ 

Proof: Applying the Arzela-Ascoli theorem several times to 



[a G C^-(^^2d) 



\a\\ck < L2, 



<L4x-y\\2yx,yeQ2d} (13) 



shows the compactness in {C^{Q,2d), \\ • Wc'')- ^'^'^ being a closed subset of (|13p proves the 
statement of this lemma. | 

Definition 5.5 (C^ - Convergence of Sets) 

r2(a„) r2(a) as n — )• 00 is defined by a„ — )• a in C'^(r22d) as n — )■ 00, where a,an and 
r2(a), r2(a„) for n G N are defined as in Definition 15.31 
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In order to apply regularity results of linear elasticity we have to require a sufficiently 
smooth boundary of the admissible domains. Therefore, additional boundary conditions on 
the design variables a are introduced which enable the construction of such domains. 

Definition 5.6 (C'^ -Admissible Domains) 

Let U^'^ be the set of design variables of Definition 15.31 and let Sjs : di}2d — ?• K be functions for 

multi-indices /3 with 1 < |/3| < k. Define U"^ = |a E If^ \ V^alan^a = ^13 G {1, A:} }• 

Choosing so that Q{a) has a C'^-boundary for every a G U"''^ the set O = {^{a) \ a G U^'^} 
denotes the family of so-called C'^-admissible domains. 

Lemma 5.7 U"-'^ is compact in (C'^(r22d), || • ||cfc). 

Proof: Note that U^'^ is a closed subset of IJ^'^, where U^'^ is already compact according to 
Lemma 15. 4[ | 



6 Existence of Optimal Shapes 

We consider optimal shape design problems in linear elasticity which in particular include shape 
optimization problems given by the cost functional (jlOp of our local and probabilistic model 
for LCF. The state problems are described by mixed problems of linear elasticity. We analyze 
a very general class of cost functionals which are not constraint by convexity assumptions and 
we optimize shapes within the family of C'^-admissible domains. As already announced in the 
previous section only a part of the boundary is subject of optimization. The abstract setting of 
Section [5] and Theorem 15.21 determine the structure of this section which leads to our existence 
results. Therefore, the technical steps will be addressed here in the following order: 

• Linear elasticity and regularity properties. 

• State equation and space given by the mixed problem of linear elasticity. 

• Convergence of functions with variable C'^-admissible domains. 

• Compactness shown by the Arzela-Ascoli theorem and by Schauder estimates. 

• Class of cost functionals and continuity properties. 

• Proof of the main existence result. 

• Applying the existence result to our probabilistic and deterministic models of fatigue. 

At first, we outline important results in linear elasticity which will influence our choice of 
state space and of the definition of convergence of functions with variable domains. Recall the 
mixed problem ([1]) which is called pure displacement problem if d^l = d^la and pure traction 
problem if dQ = BQn- But in the following we consider the so-called disjoint displacement- 
traction problem where and 90 at are disjoint. 

Since we assume a homogeneous temperature field the Lame coefficients A and /i are con- 
stant. Regularity results for the mixed problem depend crucially on the properties of the 
domain's boundary in which the elasticity equations are posed. The following statement of 
|lUj ensures the existence of a weak solution of the mixed problem. 

Theorem 6.1 (Existence of a Weak Solution) 

Let Q C M" be a domain and dVlo C dVl he measurable where dVlo has a positive area. Let 
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the Lame coefficients A, be positive constants and let f € [L^^^{ft)]^ , g G [L'^/^{d^N)f' where 
driN = dQ \ driD- Moreover, define on Vdn = {v ^ {H^{Q)]^ \ v = a.e. on 80.0} : 

B{u,v) = / XVuVv dx + / 2fi { 2. £ij{u)£ij{v) dx, 
Jn Jn y..^^ J 

L{v) = f ■ V dx + g ■ V dA. 

Jn JdUN 

Then, there exists a unique u & Vdn that satisfies 

B{u,v) = L{v) yv£VDN (14) 

and additionally J{u) = inf{J(f) | v G Vdn}, where J{v) = ^B{v,v) — L(v). 

The unique solution u £ Vdn of ([UD is called the weak solution of ([T]). The following 
inequality can be found in [7]: If 17 C M"^ is a domain the so-called Korn's second inequality 

c\HiH^(n)Y-i < ^J^^£%v)ije%v)ijdx^ = \\£%v)\\ho(^q) (15) 

holds for all v S Vdn- The disjoint displacement-traction problem has additional regularity 
properties if 90 and the forces / and g are sufficiently regular, confer Section 6.3 of |10| : 

Theorem 6.2 (Regularity of the Weak Solution for the Disjoint Displacement- 
Traction Problem) 

Let ncR^ be a domain with a C^-boundary, let f G [W^'P{n)]^ and let g G [W^^^/P^P{dn)]^ 
for some p > 6/5. Consider on Q a disjoint displacement-traction problem. Then, there exists 
a unique solution u G Vn of B{u,v) = L{v) for all v G Vn, where Vn,B and L are defined as 
in Theorem \6.1\ Moreover, u is an element of [W^''^ (0,)]^ . 

The key to the proof is to employ the fact that the previous problems (pure displacement, 
pure traction, disjoint displacement-traction problem) are uniformly elliptic and satisfy the 
so-called supplementary and complementing conditions. These conditions are introduced in 
[2] where Schauder estimates are also described in detail which can be applied to solutions of 
mixed problems. The Schauder estimates are an important ingredient in our proof of existence 
for optimal shapes in this section. 

In the following we consider on C^-admissible shaped f](a) G O the state problem which 
is given by the mixed problem, here in a rewritten form of ([T]): 



Ui = on dCl{a)D, i = 1,2,3, 'P(q^) 
, duj ( dui duj\ 

^'''l^^^.+^'l^''A^^.+^.)=3i ondn{a)N, ^ = 1,2,3, 



C is needed for Theorem 16.21 Therefore, set = 4 in the definition of O. 
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where dQ{a)D = dB{z,r) is the complete interior boundary of ^(q) , d^l{a)N = dil{a) \ 
dQ,{a)D the exterior boundarjH and where v is the normal of (9i77v(a). See Figure [3l too. 
We choose V{^l{a)) = [C^(J7(a))]^ as state space for V{a) so that we can additionally analyze 
stress gradients in our fatigue models. Therefore, we employ the following definition with q = 3 
for the convergence of functions with variable domains in O, also confer Section 2.5.2 in |2Uj . 

Definition 6.3 (C^- Convergence of Functions with Variable Domains) 

Recalling the sets O and O^^* of Definition [O] let pn : [Ci{n)]'^ [C^(J1^^*)]3 be the extension 
operator which can be derived from Lemma IA.4I for g E N \ {0}. For u E [C"^(r2)]'^ set 
M^^^t = p^u. For (Oi)iGN C 6, n e O and (ti/)igN with m E [C"?(^^^)]^ / E N, and u E [amf 
the expression n as / — )• oo is defined by ti^^* — )• n°^* in [Cq (0'^^*)]^. 

We want to find an optimal shape r2(a) E O which minimizes a functional of the form 
J{Q,u) = Jvoi(f^,u) + Jsuri^ju) with 

J^o\{^,u) = / J='^ol{x,u{x),Vu{x),V'^u{x))dx, Jsur{^,u) = / I'suri' , u{-) ,^u{-) u{-)) dA 

JQ Jan 
and u uniquely given by Q{a) as the solution of the state problem V{a). The uniqueness is 
realized by introducing appropriate assumptions with respect to the setting of 7^ (a) what is 
shown in the following. As an important result we use the Schauder estimates of Theorem 
9.3 in [2] and validate if the corresponding assumptions are satisfied. At first, we present two 
lemmas which show the existence of sufficiently regular hemisphere transformations - confer 
Definition lA.ll - and the validity of a certain inequality, respectively. These two lemmas will 
be used in the proof for the next theorem. In the following statements Banach spaces C'^''^ for 
(j) E (0, 1) occur whose definition can be found in Section 1 of [3] and in Section 7 of [2]. 

Lemma 6.4 Each $7 E O satisfies a hemisphere property where the corresponding hemisphere 
transformations are of class C^'*^ for (p E (0, 1) and have a uniform bound k with respect to O. 

Proof: Regarding Definition lA.ll we have to show that every x £ il. within a certain distance 
d > of has a neighborhood Ux with B{x, d/2) C Ux and 

UxnU = r,(SR(,)), < R{x) < 1, Ux ndn = r, (f^(,)) , 

for some hemisphere!^ '^r(x) ^-^id transformations Tx, of class C^''^. At first we consider 
x € Q within a sufficiently small distance d > of r(a) where T{a) denotes the portion of 
dQ = dO,{a) which is determined by the design variable a, see Figure HI 

Because of the definition of the basic design and of the admissible shapes a E U"''^ there 
is a C^-extension a*^^* : Q,2d^ — M of a with C ^2^* which describes a portion V^^^^a) 
of the boundary dQ,{a) beyond r(Q) and where figS* image of a C^-diffeomorphism 

T2d '■ -6(0, i?) C — )• Q,2d^ for some R > 0. This extension is needed in order to consider 
all X E r2(a) within a sufficiently small distance d > oi T(a). Now, we are able to define a 
hemisphere transformation with the required properties: 

/ {T2d{xi,X2))l \ 

r,:S^^M2^ Tx{xuX2,X-s)= i {f2d{xi,X2))2 , (16) 

V a'"*((T2d(xi,X2))l,(f2d(xi,X2))2)-X3 / 

^Note that this decomposition of the boundary depends continuously on a £ U"''^ and the two-dimensional 
Lebesgue measure of dQ{a)D is greater than a positive constant for all a € U'^'^ . 
Fr(x) denotes the flat boundary of the hemisphere E^jj;). 
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see Figure [H too. The neighborhood Ux can be chosen so that 11x^1^ = Tx{Tijii^x)) and 
Ux n do, = Tx This can be achieved by sufficiently expanding Tx{T,ji(^x)) beyond a*^^*. 

Because of the definition of the basic design and of the design variables a one can find a bound 
for the norms of the hemisphere transformations which is valid for all a G U"''^. Analogously 
the remaining hemisphere transformations can be constructed which are of a finite number and 
all have a uniform bound denoted by k. | 



The following lemma contains an inequality which can be found in Section 7 of |T]. Using 
only a few additional technical arguments a statement about the inequality's constant C can 
be added regarding its dependency on cone properties of the underlying domain 0. 



Lemma 6.5 Let Ai be a set of bounded domains in M" with a uniform cone property and let 
£ M. Then, for every e > there is a C{e) > uniform with respect to M such that 
ll''^llc"(n) ^ £||^llci(o) + ^ I'vldx holds for all v G C^{Q). 

Before applying the Schauder estimates of Theorem 9.3 in [2] to the solutions of V{q) we 
introduce the following technical notations: For arbitrary multi-indices q and (3 define 5gp to 
be one if ^ = /3 and to be zero in any other case. A vector and a matrix of multi-indices is 
given by Q{i) and I3{ij) for i = 1, . . . ,n and j = 1, . . . ,m, respectively, where each Q{i) and 
fi{ij) represents a multi-index. The A:-th component is given by Q{i)k and f3{ij)k, respectively. 

Theorem 6.6 Let the boundary value problem V{a) be given on a domain il.{a) G O. Suppose 
that the Lame coefficients are constants. Moreover, let fi G C^''^{Q'^^^) and gi G C'^''^{dQ{a)) 
for i = 1,2,3 and (p G (0,1) and let their corresponding norms be bounded by a constant kfg 
which is independent of a. Then, there are hemisphere transformations of class C^'"^ with 
uniform bound k and a unique solution u G V{Q.{a)) of the boundary value problem V{a) 
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which also belongs to (r2(a)) and satisfies 

/ 3 3 3 \ 

ll^llca.-^in) < C ( ^ ||/i||ci>'^'(c) +'^\\9h\\c^.4>{dn) + ^\Wk\\co{Q) = 1'2,3, (17) 

\i=l h=l k=l / 

for some constant C. The terms \\uk\\co(n) ca7^ be replaced by \uk\dx for k = 1,2,3 and C 
can be chosen uniformly with respect to O. 



Proof: The existence of a unique solution u € ^(^(q)) is a consequence of Theorem | 
Because Q.{a) has C^-boundary we have a unique solution u E [Vl^^'^(r2(a))]^ for arbitrary 
P > 6/5. Then, the general Sobolev inequalities - confer Section 5.6 in [T3] - will lead to 
u e [C^''^{Q{a))]^ for some (j) G (0, 1) if p is sufficiently larg^lll- 

Now, we show that the assumptions for the Schauder estimates of Theorem 9.3 in [2] 
are satisfied. Main assumptions are complementing and supplementary boundary conditions, 
uniform ellipticity with the corresponding constant A, a positive minor constant Ag^ and the 
existence of hemisphere transformation^^ of class C^''^' with corresponding norms uniformly 
bounded by a constant k, confer Sections 1, 2, 7 and 9 of [2] for a detailed description. We 
start by rewriting the equations of V{a) in the form of Yl^=i hji^, V)uj{x) = fi{x) for x G O 
and of Yl^j=i^hj{x,'^)uj{x) = gh{x) for x £ d^ljy, see (1.1) and (2.1) in [2]. Therefore, we 
write 



2 

lij{x,^) = ^ ^ aij^gix)! 



E 

\(>\=2 



(18) 
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where ■j{ij) = {5ii + 5ij,52i + S2j,d3i + 5sj) and where X]/3(j)i>o ^^^^ multi-indices 

with P{i)i > 0. Corresponding to the Neumann boundary conditions on OQn we write 

1 

Bhj{x, S) = ^ bhj,e^^ = ^ {>'^h{x)5gp(^j) + 6hjiivg{x) + ^iVj{x)5gp(^h))'E^ , (19) 
lel=o \e\=i 

where /3(1) = (1, 0, 0), /3(2) = (0, 1, 0), /3(3) = (0, 0, 1), z.(i,o,o) = ^^1,^^(0,1,0) = ^^2,^^(0,0,1) = ^^3- 
Regarding Theorem 9.3 in [2] and equations (fT8|) and (fT9|) we have Sj = O^^i = 2 and r/j = — 1 
for the Neumann condition and r^ = —2 for the the Dirichlet conditioio which implies /q = 
max{0, r/i} = and allows to choose 1 = / > /q for i, j, /i G {1, 2, 3}. 

Consider that the complementing conditions and the supplementary boundary conditions 
are satisfied, confer Section 6.3 of [lOj . The existence of appropriate hemisphere transforma- 
tions of class C^'"^ is shown in Lemma [6 . 41 where the constants d and k can be chosen uniformly 
with respect to O. According to Section 1.3 of [37] the minor constant A^q is positive and 
determined by the Lame coefficients which are constant in our case. 



'^^Note that fi and Qi, i = 1,2, 3, are continuously differentiable. 
'^^Recall Definition I A . II where also constant d is defined. 

^^Note that homogeneous Dirichlet conditions on dQ.D lead to Bhj{x,3) = S^j, gh{x) = for all x £ OQd 
and to rh = —2. 
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Now, we analyze the effect of the hemisphere transformation Tx on the ellipticity con- 
stant and follow Section 9 of [2]. Let c(V) = I]^ c^V* = I]ii,i2,... Cii.iz,...^^ • • • be 

an arbitrary linear combination of differentiation operators. It is transformed into c(V) = 
Eii,i2,...Cn,i2,...^^..., where J- = ffj^ and each Ci,,i^,,„ is a linear combination 

of the Cj^^i2,.-- with coefficients that are products of the Correspondingly, we obtain 

c(^) = c(|) with ^ = ff^Ci- According to Section 1 in [2] uniform ellipticity is described 

by the inequality < \L{x,E)\ < for ah E e M"+i and all x e Ti, where 

L{x, H) is the characteristic determinant of the PDE-system. The determinant is invariant un- 
der the hemisphere transformation in the sense that L(x,^) = L{x,S^). As the first derivatives 
of Tx and T^T^ exist and as these maps are in x uniformly bounded with respect to the norm 
II • 11(^3, <#. there is a constant oj such that a;~-'-||^||2 < HCIb < "^llClb for every x G A, ^ G ^r{x) 
and ^ = Tx{£,). Confer Chapter 1.6.2 of [18], too. Finally, this results in the uniform ellipticity 
of the transformed system with the new ellipticity constant AuP''^ . 

Applying the Schauder estimates to both cases of Neumann boundary and homogeneous 
Dirichlet conditions yields the inequality statement. A uniform choice of the constant C in 
()17p with respect to O is justified by the previous analysis of the constants d, k, Agn and A. 
The replacement of ||'Ufc||cO{o) by /q \uk\dx for k = 1,2,3 is a consequence of a uniform cone 
property of O and of Lemma 16.51 | 

By means of the following theorem and the properties of C'^-admissible domains we can 
show compactness oiQ = {(r2,n(i7)) | ^2 € O} where u{VC) uniquely solves 'P(f^). 

Theorem 6.7 There is a positive constant C such that the solution u G V{Vl{a)) of the 
previous theorem satisfies \\uj\\Q3.4,(Q-) < C for j = 1,2,3, where C can he chosen uniformly 
with respect to O and forces fi G Ci'<^(0^), G C'^^'t>{dn{a)) fori = 1,2,3, (/> G (0,1) whose 
corresponding norms are uniformly hounded hy kfg. 

Proof: The main part of the proof consists of showing that there is a constant C independent 
of G O such that ||it||_ffi(Q(a)) ^ C- We follow ideas of the proof of Lemma 2.24 in pO] : 
The weak formulation ()14p of V{a) can be rewritten in the form 

/I 3 /i 3 /I 3 

/ ^ o-ij(t;)eij(i;)(ix = / ^ /iUjdj; / '^giVidA 

where aij = Yll,i=i Cijki^ki- The ellipticity constant Ciju = Sij6kiX + fJ-{5ikSji + dudjk) is given 
by the Lame coefficients A,// and so a constant element of (7(0°^*). Moreover, the constant 
satisfies the symmetries Cijki = Cjiki = Ckiij and there exists a constant g > such that 
Cijki{x)iijiki > qiijiki for all X G Q.""^^. This results in 

Bc,{v,v)= / y2 (^ij{v)£ij{v)dx>q V eij-(7;)eij(?;)dx = g||e(?;)||Lom(a))i3 
Jn(a) f-t^ Jnia) 

for all V G Vbn- Because of the assumptions for / and g 

„ 3 ^3 



/ y^fiVidx^ / 'Y^giVidA<C\u\u^(p_(^^y^f 
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holds and the constant is independent of a S U . The previous inequahties and the weak 
equation ([UD lead to '7lk('w)ll^j:/0(Q(Q))]3 < C'||n||[j|^i(Q(ct)]S'- This and Korn's second inequality 
([15]) imply 9lk('«)||fj/o(Q(c,))]3 < C'lhll[Hi(n(a))]3 < C'l|e(^)ll[H0(n(a))]3, where C also depends on 
the constant of Korn's second inequality. As j29j shows, this constant is uniform with respect 
to a class of domains possessing a uniform cone property. Applying once more Korn's second 
inequality one obtains ||^t||/i'i(n(a)) ^ C with C independent of a. 

This result, the inequality < Y^vol(il) ||f ||_l2(q) for v € and the results of 

the previous theorem show the statement of this corollary. | 

Lemma 6.8 Let the setting of the previous theorem be given on an arbitrary sequence of 
domains (17(a„))„eN C O. For (f> G (0, 1) and i = 1, 2, 3 let fi e C^''^(?7^) and gt G C'^'^{n^) 
such that gi\dn(a)M ^ C'^''^ {d^{a) for a G U"''^ and such that their corresponding norms 
are uniformly bounded by a constant kfg. Let (a„,u„)nGN be a sequence of admissible shapes 
an G U"''^ and of the corresponding solutions Un G V{^{an)) ofV{an). Then, there exists a 

subsequence (a„j., u„^J„^^gN such that 0(a„j,) — > Q{a) and Un^ u as ^ oo for some 
a G C/"*^ and for the corresponding solution u G V{Q.{a)) ofV{a). 

Proof: Due to Lemma [STTl there is a subsequence (a^JniSN with Q{ani) — > ^{a) as — )• oo 
for some a G U"''^. According to Theorem 16.61 u^, G [C'^''^(r2(a„.j))]^ holds for every n; G N. 
Moreover, according to Theorem 16.71 there is a constant C > independent of every G O 
such that ||iinJ|[c3:'*(Q(o„;))]3 ^ Because of Lemma [A. 41 there is a constant C such that 
<f = Pn{a„^)Uni satisfies ||<f ||[p3,,A(oext)]3 < C'lhnJ|[c3,0(n{a„,))]3- Chapters L6.2 and L6.9 
of [18] show that the constant can be chosen to be independent of the design variables a. 
The previous inequalities result in a uniform bound for all ||^i^^*||j^3,,^^f^cxt)]3- unit 
ball in C^''^{Q,^^^) is compact in C^(f]°^*) due to the Arzela-Ascoli theorem we find a further 
subsequence (ar^., such that 

r2(a„J ^{a), u^^^ V as oo for some v G [C^(r2*'''*)]^. 

Because of u'^^ — > v as — s- oo the function v\ ^^^-^ solves 'P(a). According to Theorem 16.21 
V{a) has a unique solution u G V{^{a)) and v is an extension u*^^* of u. | 

As already mentioned we consider functionals which are volume or surface integrals with 
continuous integrands. In order to show continuity of the functionals we apply Lebesgue's 
dominated convergence theorem. 

Lemma 6.9 Let J-yoi, J-sur G C'^(M'^^,M) and let the set O only consist of -admissible shapes. 
For G O and u G [C^(r2)]^ with g > 1 and (j) G (0, 1) consider the volume integral 

Jyoi{^,u) = / J='yoi{x,u{x),Vu{x),V'^u{x))dx. 
Jn 

Let {^ln)n£N C O with Q,n — > ^ as n ^ oo and let {un)n£N be a sequence with Un G [(7^(0^)]^ 
and Un ^ u as n ^ oo for some u G [C'^(r2)]'^. Then, Jyoi{^n,Un) — ^ Jvoli^iU) as n ^ oo. 
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If the set O only consists of C^- admissible shapes and if we consider a surface integral 



J,ur{^,u)= / Tsur{;u{-),Vu{-),V'u{-)))dA, 

Jan 

one obtains Jguri^mUn) — )• Jsur{^,u) as n ^ oo as well. 
Proof: 

At first, we consider the volume integral. Using the characteristic function one obtains 
J(a„,u„)= / xn„ix) ■ J^voiix,u''^\x),Vu''^\x),V^u''^\x))dx. 

Because of J-Voi G C''(M^'',M) and u„ w as n — )• oo there is a constant C > such that 
the inequality \xn„{x) ■ J^miix,uT^ix),'^u'^^ix),'V'^u'^^{x))\ < C holds for ah ?i G N almost 

everywhere in Q!^^^. Moreover, i7„ — > and — t- in [Cq (0'^^')]^ as n — t- oo ensure the 
existence of 

lim (xc„(x) • J-voi(x,nr(x),V<-*(x)),V2nr(:r))) 

n— >oo 

= Xn{x) ■ J^^,i{x,u^^\x),Vu''^\x),V^u^-\x)) 

for all x £ fi^'^*. Lebesgue's dominated convergence theorem can now be applied to permute 
integral and limit: 

lim J(J^„,'u„) = lim / xnAx) ■ J'voi{x,u^n\x),Vu^^\x),V^u^^\x))dx 

n~>co n-i-cxD J^cxt 

= / hm (xnjx). J-™i(x,^zr(x),Vnr(^)),V2<-*(x))) dx 
= / J^voi{x,u{x),Vu{x),V'^u{x)) dx = J(i},u). 

With respect to the surface integral similar arguments can be used and so we only address 
technical steps which are special to integrating over a surface. Because of their definition the 
boundary of every € O is a differentiable submanifold and the surface integral 

J(0„,n„)= / Fsnr{-,Un{-),VUn{-),'^'^Un{-))dA 

J afi„ 

is the sum of integrals over a finite number |/| of map areas This number is the same 

for every G O due to the definition of the basic design and of the design variables a G U"''^. 
The integrals have the form 



-T^sur {(Pnis),Un{(Pn{s)), VUn {(Pn{s)) ,V'^ Un{iPnis))) \/ 9'^" is)ds 

where Lp\ : Jd"^ — )• is a chart and g'^"{s) the corresponding Gram determinant. The map 
93" : ^2d — ^ 1^^; ¥'n('5i) •S2) = (si , S2 ; On^si j •82))^ is the chart for the portion of the boundary 
d^n which is described by the design variablqlj an £ C^{^2d)- Because of J^sur G 
and u„ ~^ u as n — )• 00 and because of the fact that g'^" is bounded we can apply Lebsgue's 
dominated convergence theorem again to prove J{^ln,Un) — )■ J{^l,u) as — )• oo. | 



*This is the only part of the boundary which varies over the different domains of C 
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Now, we can prove the existence of optimal shapes and present our main result. 

Theorem 6.10 (Existence of Optimal Shapes for SO Problems in Linear Elasticity) 

Let F^oh^suT G C°(M''°,M) and recall 6, O of the previous definition^^. Consider the state 
problem V{a) and the functional J(i},u) = J^oi{Vt,u) + Jsur{^,u) with 



Jyol{^,u) = / J='yol{x,u{x),Vu{x),V u{x))dx, Jsur{^,u) = / J^sur{- , u{-) ,'^u{-) ,V u{-)) dA 
Jn JdQ. 

and u = u{rt) unique solution ofV{a). With respect to V{a) let the Lame coefficients X, fi 
be constant and for (p € (0,1) and i = 1,2,3 let fi £ Ci'*(0«^) and gi G C2''?^(Q^) such 
that gi\dn{a)M ^ C'^''^ {d^{o-) n) for every a € C/'"^. Moreover, let their corresponding norms be 
uniformly bounded by a constant kjg independent of a £ U"''^. Then, there exists an optimal 
shape n* €0 such that J(il*,n(0*)) < J{n,u{n)) for all G C 

Proof: Let (a„,u„)„gN be a minimizing sequence of mi{J{il.,u{Q) \Q G O)}, where Un = 
u{Qn) = u{Q{an)) is the unique solution of the state problem 'P{an)- Then, because of 

Lemma [621 there exists a subsequence {anf,,Un^)n^^^ such that f](a„j.) — > 0(a) and Un^ ~^ u 
as nfc — )• oo for some a E C/""^ and for the corresponding solution u E V{Q.{a)) of 'P(a). Since 
all prerequisites of Lemma [6.91 are satisfied we obtain J{^ri^,Unj,) — J{^l,u) as — )• oo. 
Because (a^^, , tt„^)nj.gN is also a minimizing sequence of inf{ J(r2, n(il) | E O)} the admissible 
shape r?* = r2(a) is an optimal shape. | 



Next, we apply the result of the previous theorem to the cost functional (jlOp . where now 
the number of cycles n is fixed but where the domains O = i}(a) E O are variable. 

Theorem 6.11 (Optimal Reliability) 

Let the admissible shapes ft = Q{a) of O be subject to surface driven LCF failure mechanism. 
Regarding V{a) let the assumptions of Theorem \6.10\ be imposed and let u be the unique solution 
ofV{a). Let Nfi^t be the solution of the CMB equation and consider the cost functional 



Jsur{^,u)= I — 1=—-^] dA 
an \Ndet{ea{yu{x)))' 



1 



Ndet{ea{yu{-))) 



(20) 

L™(an) 



where x i-> £aiyu{x)) represents x i— )• ROo SD ^ o o cr^(-) o Vu{x). Then, there exists an 
optimal shape $7* E O such that Jsuri^* ,ui^*)) ^ Jsuri^^ui^)) for all Q € O. 

Proof: As all prerequistes of Theorem 16.101 are satisfied we directly obtain the existence of an 
optimal shape J7* E O. | 

Theorem 16.111 proves that there is an optimal design fi* E O which minimizes the expected 
number of LCF crack initiations divided by n™. This results in a minimum for the proba- 
bility for LCF crack initiation and ensures optimal reliability. According to the remarks to 
Definition 14.11 volume driven HCF failure mechanism can be considered by volume integrals 
whose integrands only differ from the LCF integrands by the material parameters. In that 
case the existence of an optimal reliable design can be proven by means of Theorem 16.101 as 



consists of C*-admissible shapes. 
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well. Moreover, note that a possible extension of our models to an inclusion of V^n and stress 
gradients can be considered via Theorem 16. 101 

At the end of this section we prove the existence of an optimal design in case of the 
deterministic LCF model according to Definition 14.21 Note that the following result can also 
be extended to deterministic HCF. 

Theorem 6.12 (Deterministic LCF and Shape Optimization) 

Under the setting of Theorem \6.11\ consider the functional 

Jsur{^,u) = sup inf 7Vde4(ea(Vn(rj)(x))). 

QgO xedn 

There exists an optimal shape Q* ^ O such that Jsuri^* ^ui^*)) ^ Jsuri^Tu{Q)) for all G O. 

Proof: The proof of Theorem 16.101 yields all arguments needed except for the continuity 
property of Jsur- But this follows from uniform convergence required in Definition 16.31 | 

We have proven the existence of optimal designs for a general class of cost functionals 
without convexity constraints. The cases of our local and probabilistic model and of the 
deterministic model for LCF are included. Moreover, we have shown that these results can 
be extended to the case of volume driven HCF. In jl9j we address sensitivity analysis and the 
existence of shape gradients in conjunction with our probabilistic model for fatigue. 

A Appendix 

Definition A.l (Hemisphere Property, Hemisphere Transformation) 

Let A C ft he a subdomain such that dA H dfl is in the interior of T in the n-dimensional 
sense and let d > be a positive constant. If every x £ A within a distance d of dil. has a 
neighborhood Ux with Ux n dQ C F, B{x, d/2) C Ux and 

U^. n n = Tx {^R(^x) ), Uxr\dfi = Tx {Fr^x) ) , Q<R{x)<i 

for some hemisphere '^r(^x) functions Tx, T~^ of some class C^'*^, A is said to satisfy a 
hemisphere property. Moreover, the functions Tx are called hemisphere transformations. 

Definition A. 2 (C^ ^f'- Boundary, Lipschitz Boundary) fl8^ 

Let C be a bounded domain, let fc G N and < (f) < 1. The subset J7 has a C'^''^-boundary 
if at each point xq G dfl there exists B = B{xq, r) for some r > and an injective C'''''^-map 
if; : B ^ D CW such that ^(^ n 1^) C M!;: = {x G M" | x„ > 0}, tpiB n dQ) C dR^ and 
^/^~^ G C^''^{D). If the maps tp are only Lipschitz continuous Vl has a Lipschitz boundary. 

Remark A.3 

Let be a bounded domain in R"" and let A; G N and < </) < 1. If every xq G dVl has a 
neighbourhood in which the boundary is locally described by a graph of a C'^''^-function of 
n — 1 of the coordinates xi, . . . ,Xn the domain fl has a C'^''^-boundary. Note that the converse 
is true if A; > 1. 

Lemma A. 4 (Extension Lemma) f7^ . Part I, Section 6 

Let f]''^* C he open and let ft be a C''''^ -domain with Tl C rJ*"^*, k > 1 and < < 1. For 
(p = it is C^'^ = . If u £ C^''^{Q) there is a function w G CQ''^(il'^^*) such that w = u in fl 
and \\w\\(jk.,^(^Qext-^ < C'||ti||f;fc,</)(Q) for a constant C depending only on k,Q, and fl^^^. 
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